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We present preliminary results for Bk calculated using improved staggered fermions with a mixed 
action (HYP-smeared staggered valence quarks and AsqTad staggered sea quarks). We investigate 
the effect of non-degenerate quarks on Bk and attempt to estimate the ff(a 2 ) effect due to non- 
Goldstone pions in loops. We fit the data to continuum partially quenched chiral perturbation 
theory. We find that the quality of fit for Bk improves if we include non-degenerate quark mass 
combinations. We also observe, however, that the fitting curve deviates from the data points in 
the light quark mass region. This may indicate the need to include taste-breaking in pion loops. 



XXIV International Symposium on Lattice Field Theory 
July 23-28 2006 
Tucson Arizona, US 



* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http://pos.sissa.it/ 



Bk in unquenched QCD 



Jongjeong Kim 



1. Introduction 

Indirect CP violation in the neutral kaon system has long been established experimentally. It 
is parameterized by e, which is measured to high precision. In the standard model, e is given by 



e =C e exp(/7r/4) ImA, X B K (1.1) 

X = RsX c [t]iSo{x c ) - mS 3 (x c ,x t )] - ReA,r/ 2 5 (x ; ) (1.2) 

h = V*V id , x i = m 2 i /Ml (1.3) 

~ E 6V2n 2 AM K 



where is the CKM matrix element and Sj are the Inami-Nam functions. The kaon bag parameter, 
Bk, is defined by 

K l(Ko\sr^r 5 d\o)(o\sy^r5d\Ko) 

B K = C{pL)B K {pL) (1.6) 

TO 

C(ji) = a s (n) 2 "o [l + a s (n)J 3 ] (1.7) 

where B% is renormalization-group invariant. Hence, a precise determination of Bk constrains the 
CKM matrix — any violation of the unitarity of this matrix would require new physics beyond the 
standard model. Therefore, there has been a long-standing attempt to calculate Bk using lattice 
QCD and other methods. 

Staggered fermions have two advantages for a calculation of Bk- First, they preserve an axial 
U(l) symmetry which constrains the desired matrix element to have a chiral expansion similar 
to that of the continuum, and, second, they are cheap to simulate compared to formulations with 
exact (or almost exact) chiral symmetry (domain wall and overlap fermions). A drawback is the 
need to take the fourth-root of the fermion determinant when generating configurations, leading to 
unphysical effects for a / 0, and the need to use a fitting function containing such effects. Clearly 
one must also assume the validity of the rooting trick in the continuum limit. 

Unimproved staggered fermions, however, are not suitable for a precision calculation. They 
suffer from (1) large scaling violations, (2) large taste symmetry breaking, and (3) large perturbative 
corrections in matching factors. In order to alleviate these problems, there have been a number of 
proposals to improve staggered fermions. Two of them have been used extensively in the lattice 
community: (1) AsqTad staggered fermions and (2) HYP/Fat7 staggered fermions. It turns out that 
HYP staggered fermions reduce the taste symmetry breaking [jl|] as well as one-loop perturbative 
corrections [^] more efficiently than AsqTad staggered fermions, and have much smaller scaling 
violations [^]. In other words, HYP staggered fermions offer significant advantages over AsqTad 
staggered fermions. 

Here we perform a numerical study using valence HYP staggered fermions. We use the un- 
quenched MILC lattices (Np =2 + 1), in which the sea quarks are AsqTad staggered fermions, so 
that we are using a "mixed action". The details of the simulation parameters are given in Table [[]. 

Although the one-loop contributions to matching factors with HYP fermions are small, the 
unknown two- and higher-loop contributions lead to a significant uncertainty when one aims at a 
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parameter 


value 


P 


6.76 (unquenched QCD) 


sea quarks 


2+1 flavor AsqTad staggered fermions 


valence quarks 


HYP staggered fermions 


l/a 


1.588(19) GeV 


geometry 


20 3 x 64 


# of confs 


640 


sea quark mass 


am u = anid = 0.01, am s = 0.05 


valence quark mass 


0.01, 0.015, 0.02, 0.025, 0.03, 0.035, 0.04, 0.045, 0.05 


7 


ps 1.0 



Table 1: Parameters for the numerical study 



precision calculation. We aim to address this problem using both two-loop calculations and non- 
perturbative renormalization. 



2. B K for degenerate quarks (m x = m y ) 

This numerical study is a follow-up to the work in Ref. [01]. We place U(l) noise sources at 
t = and t = 26. These couple to the pseudo-Goldstone pion (spin-taste 75 <g) (§5), while projecting 
against all non-Goldstone pions. When we use this U(l) noise source (at t = 0) for two-point 
correlation functions, we observe a noticeable contamination from excited states when < t < 10. 
Hence, we choose the fitting range to be 10 < t < 15 in order to exclude contamination from both 
sources. 

In this study, we fit our data to the prediction of continuum partially quenched chiral perturba- 
tion theory. The result for Nf = 2 + 1 flavors was given in Ref. Ph. For degenerate valence quarks 
(m x = my) it takes the form 



ci 1 + 



1 



487T 2 / 2 

.2 ' 



^conn ^disc 



6<r( m 2)-12/(m 2 ), 



+ C2m xy + c 4 mi y , 



^disc 



0, 



(2.1) 
(2.2) 



where /= 132 MeV and 



l(X)=Xlog(^), l(X) 



up to finite volume corrections. The notation for meson masses is 



(2.3) 



m u = 2\lm u , 
m\ = 2\xm x , 



lily 



s — 2\Lm s , = (my + 2m s )/3 

2\im y , m 2 xy = ll (m x + m y ) 



(2.4) 
(2.5) 



where m u = 7^ m s are sea quark masses and m x , m y are valence quark masses. 
In Fig. [IJ we plot Bk as a function of M K for degenerate valence quarks (m x 



m y ). We fit 



the data to the form in eq. [2.1[ with the cut-off scale chosen to be A = 4nf. The fitting results 
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Figure 1: Bk versus m\ for degenerate quark mass combinations (m x = m y ) 

are summarized in Table ||. The X 2 is rather high, considering that the data are correlated and we 
do not include the full covariance matrix. In particular, we observe that the curve does not fit the 
lightest data point very well. 

One possible explanation for the mismatch as small quark mass is that it is due to taste- 
breaking in pion loops. As noted in Ref. [§], non-Goldstone pions dominate the loop contributions 
which lead to curvature at small m 2 K , and this reduces the expected curvature away from the contin- 
uum limit. This prediction is qualitatively consistent with what we observe. However, this behavior 
could also due to the finiteness of the volume, which also becomes important only for small m 2 K . 
Clearly this needs further investigation. 



3. Bk for non-degenerate quarks (m x ^ m y ) 

A more stringent test of the applicability of chiral perturbation theory is provided using non- 
degenerate quarks. The continuum prediction is [9] 



B K = cA\ + 



1 



^conn ~i~ ^dis 



487T 2 / 2 

The contribution from quark-connected diagrams is 



(mjf — my 



2\2 



ill 



AT 



6m%l{m%) - 3/(m|) (l + - 3l(m 2 Y ) (l + ) , 



(3.1) 



(3.2) 
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parameters 


unit 


average 


error 


C\ 


1 


0.4422 


0.0145 


C2 


GeV- 2 


-1.4616 


0.1611 


C3 


GeV- 2 






C 4 


GeV- 4 


1.4674 


0.1913 


£ 2 /d.o.f. 


1 


0.5396 


0.3 



Table 2: Fitting results of Bk for degenerate quark mass combinations (m x = m y ) 



while that from diagrams involving a hairpin vertex is 

/disc = (Ix+h+I^/ml (3.3) 



-xy 

57 i,( m % + m x){ m l- m x){ m l- m x) w 1 Am\^ + m\)(m 1 1 -m\)(m\-m\ 
K m x) : 2 2 l ( m X) " r~2 2\2 



mt — m£ ( m„ — m% 



-l(m\) 



2(m%+m\)(ml-m\)(ml-m\) 



' in 



2 -m\)(m\-m\) 



-/(mi-) : - 2 : ( 3 - 4 ) 



2 m|) (m| - m|) - {ml. + m|)(m| - m|) - (mL + m x ){ml - 



/y=/x(m|^m|) (3.5) 
; _ i ( m 2 ^ ( m I ~ m r ) 2 ("4 + m \ ) ~ OT g ) ( m s ~ m ^ ) (36) 

In Fig. ||, we plot 5^ data as a function of M\, including data for both degenerate and non- 



degenerate valence quarks. We fit the data to the form in eq. 3.1. The results of the fit are summa- 
rized in Table ||. 

When we compare Table ^ with Table we observe the following: (1) # 2 /d.o.f. is much 
smaller when we include non-degenerate combinations; (2) the low energy constants c\, C2 and C4 
are consistent within statistical uncertainty between the two fits; and (3) for our choice of the cut- 
off scale, A = 4nf, the low energy constant C3 is negligibly small. This comparison tells us that, 
for this cut-off, the effect of non-degenerate quark masses are dominated by the chiral logarithmic 
terms contained in 7 conn and /disc, whereas the C3 term plays no role. Overall, we learn that, while 
we can in principle determine c\, C2 and C4 from the degenerate combinations alone, we have more 
confidence in the result when we include non-degenerate combinations. Furthermore, if we are to 
include higher order terms in chiral perturbation theory, as might be required at the heavier quark 
masses, the additional non-degenerate points will be essential. 

The primary reason to study non-degenerate valence quarks is in order to extrapolate to the 
physical kaon, in which the non-degeneracy is almost maximal. The data shows a systematic 
increase in Bk with increasing non-degeneracy. This effect is statistically significant though nu- 
merically small. When we have a satisfactory description of our data we will be able to extrapolate 
to the physical kaon mass. 
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Figure 2: Bk versus m 2 K for both degenerate and non-degenerate quark mass combinations. The fit line is 
for the degenerate combinations only. 

4. Staggered chiral perturbation theory 

A shortcoming with our present analysis is that it assumes that Goldstone and non-Goldstone 
pions are degenerate. We know, however, that taste symmetry is broken at &(a 2 ), leading to split- 
tings between pions with different tastes. A systematic way to incorporate such ff(a 2 ) effects is 
provided by staggered chiral perturbation theory [§,0]. In the case of Bk, Van de Water and Sharpe 
have determined the next-to-leading-order (NLO) prediction of staggered chiral perturbation the- 
ory l@l ■ ( m f act > a small modification of this result is necessary to describe our data, since we 
are using a mixed action.) Unfortunately, the result contains 21 unknown low energy constants, 
even at a single lattice spacing, and these must, at present, be determined from a fit to the lattice 
data. (Some can be determined in the future by calculating other matrix elements, or using non- 
perturbative renormalization to determine matching factors.) This is why we have chosen so many 
mass combinations (45 at present). 

The fitting is further constrained if we know the masses of the non-Goldstone pions. These 
enter into loops, and, as noted above, lead to a modification of the curvature at small m 2 K . Because 
of this, in a companion work, we are calculating all pion masses (as well as those of other hadrons) 
using HYP-smeared valence staggered fermions [j|]. A striking result of that study is that the 
taste-breaking effects are very small, much smaller than those with AsqTad improved staggered 
fermions. Thus it may be that we can get away with a continuum-like fit except at the smallest 
quark masses. 

An exception to this comment concerns the contributions from quark-disconnected diagrams 
involving hairpin vertices. These bring in taste-singlet pions composed of (AsqTad) sea-quarks, for 
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parameters 


unit 


average 


error 


c\ 


1 


0.4451 


0.0149 


Cl 


GeV- 2 


-1.5049 


0.1672 


C3 


GeV- 2 


0.0305 


0.0083 


c 4 


GeV- 4 


1.5366 


0.2021 


X 2 /d.o.f. 


1 


0.1888 


0.0935 



Table 3: Fitting results of Bk for non-degenerate quark mass combinations (m x ^ m y ) 

which (a 2 ) terms are known to be large. Fortunately, the masses of these pions are available from 
the MILC collaboration studies. 

As always when doing chiral extrapolations, it is important to push to as light a quark mass 
as possible (within the constraints of finite volume effects and computer time). In particular, this 
allows one to more thoroughly test the applicability of chiral perturbation theory at the order we 
work. Because of this, and to provide more data for fitting, we are presently extending the calcula- 
tion to a lighter valence quark mass (m = 0.005), which will also increase the total number of mass 
combinations to 55. 
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